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ABSTRACT 

Quasi-periodic propagating intensity disturbances have been observed in large coro¬ 
nal loops in EUV images over a decade, and are widely accepted to be slow magne- 
tosonic waves. However, spectroscopic observations from Hinode/EIS revealed their 
association with persistent coronal upflows, making this interpretation debatable. We 
perform a 2.5D magnetohydrodynamic simulation to imitate the chromospheric evapo¬ 
ration and the following reflected patterns in a flare loop. Our model encompasses the 
corona, transition region, and chromosphere. We demonstrate that the quasi periodic 
propagating intensity variations captured by the synthesized Solar Dynamics Observa¬ 
tory / Atmospheric Imaging Assembly (AIA) 131, 94 A emission images match the pre¬ 
vious observations well. With particle tracers in the simulation, we conhrm that these 
quasi periodic propagating intensity variations consist of reflected slow mode waves and 
mass flows with an average speed of 310 km/s in an 80 Mm length loop with an average 
temperature of 9 MK. With the synthesized Doppler shift velocity and intensity maps 
of the Solar and Heliospheric Observatory/Soiax Ultraviolet Measurement of Emitted 
Radiation (SUMER) Fe XIX line emission, we conhrm that these rehected slow mode 
waves are propagating waves. 

Subject headings: magnetohydrodynamics(MHD) — Sun: corona — Sun: hares — Sun: 
oscillations 


1. Introduction 


The study of magnetohydrodynamic (MHD) waves in the solar atmosphere is an indepen 
dent tool to understand the energy release processes, particle acceleration or heating mecha¬ 


nisms and to diagnose the pla, s ma parameters indi rectly by coronal seismology (IRoberts 


200C 


De Moor tel &: Nakariak(^ 2012 : Liu fc Ofman 2014). MHD seismol ogy was successfully applied 
i n estimating the coronal magnetic held dNakariakov fc Ofman 2001 ). transverse loop structuring 


(Goossens et al 


2002 


(Van Doorsselaere et al 


Aschwanden et al.ll2003 ), polytropic index and ther mal conduction coe fficient 
201 Ibl ). and tlie magnetic topology of sunspots ( Yuan et al. 2014al E]). 


Standing longitudinal slow-mode oscillations were hrst discovered in the Doppler shift of hot 
emission lines (i.e., Fe XIX and Fe XXI) with formation temperature greater than 6 MK, by the 
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Solar Ultraviolet Measurements of Emi tted Radiation (SUMER) spectrograph onboard the Solar 
and Heliospheric Observatory (SOHO) ( Wang et al. 2002 . 2003a b; Wang 2011 '). Similar Doppler- 
shift oscillations have been detected by Yohkoh /BCS in even hotter emission lines of S XV and Ca 
XIX, with formation temperature 12~14 MK (jMariskal l2005l . l2006l ). The oscillations are strongly 
damped within a couple of periods and are nsually observ ed in association with the soft X-ray 
brightenings or even up to M-class flares (jWang et al.l 120071 ). 


Excitation of slow magnetoacoustic oscillations in hot coronal loops has been intensively stndied 
theoretically. The compressible natnre of the longitudinal oscillations and their long periods led to 
their interpr etation in terms of sta nding slow magnetoaconstic oscillations damped dne to thermal 


conduction (jOfman &: Wangl l2002l i . In order to explain the observed damping time of the oscilla- 


effects (Nakariakov et al. 2004: 


Tarovan et al 


20071 : 


Tsiklanri et al 


Gruszecki &: Nakariakov 


Tarovan et al. 

2Q05; 

Selwa et al. 

2Q05i 

2007; 


201lh . accounting for viscosity, multi-dimensional ge¬ 


ometry, stratihcation, nonlinear steepening, and mode coupling. 

Quasi-periodic pulsati ons (QPP) observed in solar and stellar flares have b een intensively 
studied for several decades (jNakariakov &: Melnikov! l2009l : I Anhnogentov et al.ll2013l i. The origin of 
QPPs still remains unclear, but one of the widely accepted theories is the modulation of QPP by 
magnetohydrodynamic (MHD) oscillations. Short period (snb-minute) oscillations are believed to 
be induced by fast mode waves, whil e those with periods of tens of seconds are ascribed to modnla- 
tions by slow mode MHD waves (e.g. IVan Doorsselaere et al.ll2nilal ). Coronal MHD oscillations are 
directly seen in various bands with modern instruments with high temporal and spatial resolution, 
which provide researchers with MHD wave diagnostics to identify physical conditions in flaring sites 
and mechanisms operating in them. 

Recently, the Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA) pro¬ 
vided high temporal and spatial resolution observations of slow mode oscillations in the solar 
corona. The hrst simnltaneq n s obs ervations of the electron density and EUV intensity oscillations 
were reported bv iKim et al.l ( 20121 ) through Nobeyama 17 GHz and AIA 335A channels, respec¬ 
tively. Kumar et al. (2013) reported the first direct observation of a propagating EUV disturbance 
(i.e., slow mode wave) in hot coronal arcade loops captured only in the AIA 131 and 94 channels. 
The wave was excited by an impnlsive flare which occnrred at one of the footpoints of the arcade 
loops. It showed multiple reflections between the opposite footpoints of the arcade loops (Knmar 
et al. 2013, 2015). The observed properties of these oscillations matc h the SUMER Doppler-shift 
oscillations associated with the slow magnetoacoustic mode. However, IWand (120111 ') interprets the 
SUMER Doppler-shift oscillations as standing slow waves due to the associated intensity variations, 
which show roughly a quarter-period phase delay to the Doppler signal in some cases. 


In this paper, we investigate the slow magnetoacoustic oscillations in a flare loop by a 2.5D 
magnetohydrodynamic simulation. The paper is then organized as follows: in §2 we describe the 
numerical setup; in §3 we show resnlts of simnlations and discnss the details; and conclnsions are 
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drawn in 54. 


2. Computational Aspects 
2.1. Governing Equations and Initial Setup 


Our numerical setup includes gravity, anisotropic thermal conduction and radiative cooling 
and parametrized heating terms, in a domain of -40 Mm < x < 40 Mm and 0 < y < 50 Mm. The 
governing equations are as follows: 

dp 


dt 


+ V • (pv) = 0, 


djp^) 

dt 


-h V • ( /9VV Ptotl - 


BB 


= PS, 


dE 


-Ev -h ptot^ - 

SB 


B 


Po 


B = pg • V + V • (k ■ VT) - Q + 17, 


_ + V • (vB - Bv) = 0, 


( 1 ) 

(2) 

(3) 

(4) 


where T, p, B,v, and I are respectively temperature, density, magnetic field, velocity, and unit 
tensor. The total energy density is E = p/ (7 — 1) -|- pr'^/2 -|- B^/2po and the total pressure is 
Ptot = p + B^/2po; g = PoRq/ (Rq + y)‘^y is the solar surface gravitational acceleration with po 
as —274m/s^; H and Q are respectively the heating and radiative loss terms; k is the thermal 
conductivity tensor. Assuming a 10:1 abundance of hydrogen and helium of completely ionized 
plasma, we obtain p = l.lmpren, where is the proton mass and tt-h is the number density of 
hydrogen. We use the ideal gas law p = 2.2,n}ik-BT with the ratio of specific heats 7 = 5/3. W e 
adopt Q = 1.2n^A (T) as the radiative loss function for optically thin emission ( Colgan et al.ll2008l ). 
Below 10,000 K, we set A (T) to vanish because the plasma there is optically thick and no longer 
fully ionised. The term containing k = «:||bb quantifies the anisotropic thermal conduction along 
the magnetic held lines with the Spitzer conductivity «:|| as 10“®r®/^ erg cm“^ s“^ The 

hux of anisotropic thermal conduction has a ceiling, -sign(Vr)5(^pcf as the saturated hux, and Cg 
is the isothermal sound speed. The correction factor = 1 is set according to the values suggested 
for the coronal plasma (jGiulianilll984l . iFadevev et al.ll2002l . and references therein). 


We employ a linear force-free magnetic held for the initial magnetic conhguration, which is 
characterised by a constant cut-of-plane angle Oq as follows: 

D o f■ n ( 7rysin6»o 
Bx = -Bo cos I — I smPoexp I- - - 

R R • ('^x\ f Trysin^o 

'Kx\ „ f Try sin 00 
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with 9q = 30°, the angle between the arcade and the neutral line {x = 0,y = 0). Lq = 80 Mm is 
the horizontal size of our domain, and we adopt Bq = 50 G. 

For the initial thermal structure, we set a uniform temperature of 10,000 K below a height of 2.7 
Mm and choose a temperature profile with height ensuring a constant verti cal thermal condu ction 
flux (i.e., ndTjd y = 2 x 10® erg cm“^ s“^) above this height as used in iFang et ahl ( 20131 ) and 
Xia et al.l (|2012l b The initial density is then derived by assuming hydrostatic equilibrium with a 
number density of 1.2 x 10^® cm“^ at the bottom and the initial velocity field of all plasma is 
static. We employ a background heating rate decaying exponentially with height into the whole 
system all the time, Bq = cgexp (—^) where cq = 10“^ erg cm“^ s“^ and Aq = 50 Mm. This 


heating is meant to balance the radiative losses and heat conduction related losses of the corona 
in its equilibrium state. With the above initial setup, the whole system now is out of thermal 
equilibrium. Therefore, we integrate the governing equations until the above configuration reaches 
a quasi-equilibrium state at 72 minutes after initialisation. Then we reset the time of the system 
back to zero for the next stage of simulation. As a result, the final relaxed state of the system 
is identified as the time when the maximal residual velocity in the simulation is less than 5 km 
s“^. Panel (a), (b) and (c) in Fig. [T]show the number density, temperature and AIA 131 A of the 
relaxed system, respectively. 

We use the MPI-parallelized Adaptive Mesh Refinement Ver satile Advection Code MPI-AMRVAC 


(jKeppens et al 


2012 


Porth et al.ll2ni4l : lKeppens fc Porthll2ni4l ) to run the simulation. An effective 


resolution of 1024 x 640 or an equivalent spatial resolution of 79 km in both directions is obtained 
through four AMR levels. Considering the left and right physical boundary, density, energy, y and 
z momentum components. By and Bz are set as symmetric, while Vx and Bx are taken antisymmet¬ 
ric to ensure zero face values. In the bottom boundary ghost cells, we use the primitive variables 
(p, v,p, B) to set all velocity components antisymmetric to enforce both no-flow-through (vertical) 
and no-slip (horizontal), while the B are fixed to the initial analytic expressions of equation ([5]), 
and the stratification of density is kept at pre-determined values from the initial condition, as well 
as the pressure. For the top conditions, we set all velocity components as antisymmetric, and adopt 
a discrete pressure-density extrapolation from the top layer pressure with a maximal temperature 
Ttop = 2 X 10® K. 


2.2. Imaging and Spectroscopic Modelling 


To synthesize the observational features of SDO/AIA channels, we calculated the AIA tem¬ 
perature respon se function i Fnfnp , T)[DN cm® s“^]. The detail of the forward modelling method, 
can be found in lYuan et al.l (|2015l l. The source of the forward modelling code (FoMo) is available 
at https://wiki.esat.kuleuven.be/FoMo. Then we assume that the flux Fa{x,y)[ DN s“^] will be 
integrated along the LOS for a width of W =1 Mm, 


Fa{x,y) = Ka{ne,T)nl x W. 


( 6 ) 
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We synthesized the AIA 94, 131 A channel emission which could image the flare loop at 6.4 MK 
and 10 MK, respectively, and the AIA 304 A channel emission line which would represent the 
transition region and top chromosphere (0.05 MK). 


Furthermore, we synthesized the extreme ultraviolet (EUV) emission intensity I\g [ergs cm“^] 
of a specihc spectral line Aq for optically thi n plasma along LO S for a width of W =1 Mm. The 
detail s of the method can also be found in I Yuan et al.l (|2015l l and Antolin fc: Van Doorsselaere 
(|2013l h and the intensity is given by 


ho = ^Gxo{ne,T)nlxW, 


(7) 


where Ab is the abundance of the emitting element relative to hydrogen, and Gaq [ergs cm“^ s“^ 
is the contribution function that contains the terms relative to atomic physics, as a look-up table 
for SUMER Ee XIX 1118.1 A line in which most spectroscopic observations of standing slow waves 
were performed (|Wandl201lh . 


2.3. Triggering the Flare and Chromosphere Evaporation 


The flare is triggered by a finite duration heat pulse defined by the function Hi located at 
the right footpoint between x = 22, 24 Mm as the formula below (note that the heating rate H in 
equation (3) is H = Hq + Hi). In our simulation, the heat pulse is controlled by f{t), which starts 
at time t = 0 and switches off at t = 180 seconds. The energy input by Hi is around 3 x 10^® erg 
s“^, with an assumed thickness of 1 Mm along the third axis. This energy input is suitable for a 
normal solar flare energy release. The energy of Hi is quickly transported by thermal conduction 
to plasma at the footpoints. 


Hi = ciexp{-{y-yc)^/X^)f{t) if A{xi,0) < A{x,y) < A{x 2 , 0 ) (8) 
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( 10 ) 


where ci = 16 erg cm“^ s“^, yc = H Mm, = 10 Mm^, xi = 24 Mm, and X 2 = 22 Mm. The heat 
pulse is located close to the loop footpoint, i.e. ~ 0.3 Mm above the transition region mimicking the 
footpoint heating by dissipated non-thermal particles. A{x, y) is the magnetic potential depending 
on the location and decaying exponentially with height. Because the magnetic potential along a 
single magnetic field line is constant, we only add extra heating Hi at one feet of a magnetic flux 
tube consisting of the magnetic field lines identified by A{x,y) in the range of xi < x < X 2 . The 
large ratio between ci and cq as 1.6 x 10^ highlights the extremely violent energy release of the 
solar flare. 













3. Results and Discussion 


With the number density and temperature maps from the simulation, and the methods briefly 
presented in §2.2, we can calculate the synthesized emission maps for both AlA and SUMER. We 
now discuss this simulation in more details. The three columns in Fig. [1] display the temporal 
evolution of number density, temperature and synthesized AIA 131 A emission intensity maps 
at t « 0, 83, 166 and 581 seconds, respectively. The extra flare heating Hi at the right footpoint 
(x = 23 Mm) imitates the explosive energy release with accumulated relativistic particles, suddenly 
dissipated in the upper chromosphere and transition region. This enormous energy heats the cold 
chromospheric plasma (around 0.02 MK) to an average temperature around 10 MK, as presented 
by the panel (e). The heated plasma is strongly evaporated into the confined loop as shown in panel 
(d). The velocities of the plasma are represented by the arrows in panel (d), (g) and (j), which 
show the directions of plasma movement. Compared with panel (d) and (e), panel (f) indicates that 
the synthesized emission in AIA 131 A channel is dominated by the density distribution rather 
than the temperature distribution. The front of the evaporated hot flow violently impacts the left 
footpoint (x = —23 Mm) at t ~ 163 seconds with a speed up to 600 km/s. After the violent impact 
on the chromosphere, a reflected pattern rises up towards the loop apex as shown in panel (g), 
(h) and (i) at t ~ 166 seconds. The AIA 131 A emission is mainly sensitive to 10 MK plasma. 
Still, there are subtle effects noted from comparing panel (h) and (i): the left footpoint of the loop 
has a strong synthesized AIA 131 A emission in panel (i), but this would not be expected from 
the local temperature in panel (h). This strong emission is caused by the high compression and 
heating in the left footpoint, due to the wave impact. At t ^ 446 seconds, the reflected pattern 
spreads back to the right footpoint (x = 23 Mm). Then panel (j) at t ps 581 seconds represents the 
second reflection rising from the right footpoint and shows that the average density in the whole 
confined heated loop is clearly higher than at the beginning (panel (a)). Although the maximum of 
temperature in panel (k) decreases to around 8.5 MK, we still could observe the reflected pattern 
in AIA 131 A emission in panel (1). An animation of Fig. [T] is available. 

In order to quantify and study these reflected patterns, we extracted a loop drawn by the 
white lines in panel (f) of Fig. [TJ This fixed loop is defined by a field line identified at the start 
of the simulation, with a fixed width of 1200 km. We plot average values of the number density, 
temperature and synthesized AIA 131 A emission inside this loop to time-distance maps in panel 
(a), (b) and (c) of Fig. [2l respectively. The zero and end points along s (vertical) axis in Fig. [2]mean 
the left footpoint (x = —23 Mm) and right footpoint (x = 23 Mm) in the domain, respectively. The 
reflected patterns are clearly seen as ridges in all three panels. We add a particle tracer at t ~ 166 
seconds to trace and observe the movement of the plasma, especially when the reflected patterns 
sweep over it. The location of the particle tracer (x = 0, y = 25 Mm) is marked by a white cross in 
panel (i) of Fig. [H The temporal evolution of the particle tracer movement is shown as the dotted 
line in panel (a) of Fig. [2j Because the third row of Fig. [1] shows the moment when the front of 
the initial excited flow impacts the left footpoint (s = 2 Mm) and the particle tracer is on the tail 
of the flow, the initial velocity of the particle tracer is towards the left footpoint (s = 2 Mm). The 
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particle moves downwards in panel (a) of Fig. [2] until it is swept over by the reflected pattern at 
t ~ 235 seconds (this is tl), then it turns around towards the right footpoint (s = 70 Mm). As 
shown in panel (a), the path of the particle displays three turnings (tl, t2 and t3), meaning that 
it is swept over by the reflected patterns three times. Although the particle is clearly in sync with 
the reflected patterns, the speed of the tracer is slower than the reflected patterns. The particle 
tracer behaves like a pendulum, where it is seen to oscillate back and forth, induced by the reflected 
patterns. There is another online animation of evolution of virtual particles in panel (a) of Fig. [2l 

One important characteristic used by observational studies to identify this kind of solar prop¬ 
agating disturbances as a slow mode magnetosonic wave is the agreement between the estimated 
coronal sound speeds and speeds of the reflected propagating disturbances. Based on our simu¬ 
lation, we use a more accurate method to verify this agreement. The six red solid lines in panel 
(b) and (c) of Fig. [2] show the paths of an imaginary particle propagating with the local sound 
speed, calculated by the local plasma temperature. We could consider these lines as paths of sound 
waves. In panel (c), line 1 represents the route of a sound wave initiated from the right footpoint 
(s = 70 Mm) at the same time when the extra heating Hi starts. We find that line 1 has a perfect 
agreement with the initial excited flow from the right footprint (s = 70 Mm) to the left footpoint 
(s = 2 Mm) in panel (c). This agreement indicates that the initial excited flow propagates with 
the sound speed. However, in panel (b) we find that the temperature rises earlier than the arrival 
of the flow (line 1). That is because of the faster propagation speed of the thermal conduction 
discontinuity, which also introduces a weak evaporation to slightly increase the density at the left 
footpoint (s = 2 Mm) around f « 1 minutes as shown in panel (a). Line 2 represents the path 
of a sound wave propagating from the left footpoint (s = 2 Mm) to the right footpoint (s = 70 
Mm), assumed to set out when the excited flow propagating along line 1 impacts the left footpoint 
(s = 2 Mm). Fig. [2] show that there are AIA 131 A emission, density and temperature changes at 
the top end of line 2, the right footpoint (s = 70 Mm). These increments indicate that a reflected 
wave propagates along line 2, and impacts to trigger another reflected wave rising from the right 
footpoint (s = 70 Mm). The analysis of line 1 and line 2 confirm that these reflected patterns have 
a wave component. 

However, unlike line 1, the synthesized AIA 131 A emission of the reflected pattern from the 
bottom end of line 1 does not behave similarly with line 2 in panel (c) of Fig. [21 The reflected pattern 
propagates slower than the sound wave along line 2, especially the part from the left footpoint {s = 2 
Mm) to s = 30 Mm. Because the initially excited flow from the right footpoint (s = 70 Mm) is 
triggered by a finite duration heating pulse, it relates to flows initiated within the same time range 
of Hi from t = 0 to t = 180 seconds, which is clearly observed in all panels of Fig. [2l As a result, 
the rising reflected pattern at the left footpoint (s = 2 Mm) from t ~ 161 seconds encounters the 
rest part of the initially excited flow which still propagates towards the left footpoint (s = 2 Mm). 
Panel (c) shows that this collision delays the propagation of the reflected pattern, indicating that 
both the excited flow and the reflected pattern contain a mass flow component. Line 3 in panel (c) 
is another path of a sound wave which arrives at the right footpoint (s = 70 Mm) simultaneously 


with the main part of the reflected pattern. The middle piece of Line 3 shows that after passing 
through the “collision” region, the reflected pattern propagates with the sound speed again. As 
well as for line 1, line 3 in the panel (b) temperature map shows that the thermal conduction 
discontinuity propagates faster than the sound wave. The behaviours of lines 1, 2 and 3 in panel 
(b) and (c) indicate these reflected patterns contain both wave and mass flow component. This is 
also confirmed by the line 4, 5 and 6. All of them tell the same story that the patterns observed 
in the synthesized AIA 131 A channel emission is dominated by a repeatable wave component, 
and modulated by the mass flows where collisions can temporally retard or redirect actual mass 
flows. The wave-flow behaviour of the excited disturbance is reminiscent of shadow water waves 


(e.g. IKundu et al.l l2014l i for which a single pulse also gives a longitudinal displacement to trace 
particles. Another interesting phenomenon is that the highest temperature in our simulation is not 
produced initially by the flare heating Hi, but by the collisional compression of two patterns as 
shown in panel (b). 


The slow magnetoacoustic waves observed by iKumar et all (120151 1 are thought to be propa¬ 


gating waves, rather than standing waves as observed bv IWand (j2nii[ L The standing waves have a 


unique characteristic, a quarter-period phase delay between the associated intensity variations and 
the Doppler signal. In Fig. [3l panel (a) shows a top slit view of the Doppler shift in the synthesized 
SUMER Fe XIX line, and panel (b) shows a top slit view of the intensity of the synthesized SUMER 
Fe XIX line emission. A top slit view means visual LOS from the top of the system {y = 50 Mm) to 
the bottom of the system (y = 0 Mm). The top view of the Doppler shift in panel (a) shows clear 
reflected patterns, and so does the intensity map in panel (b). We extract the black line (x = 10 
Mm) in Fig. [3] to compare the synthesized Doppler signal and the associated intensity variations 
and identify which kind of wave we have in the simulation. Although the top slit views of the 
synthesized SUMER Fe XIX line emission integrate the quantities from the top to the bottom, the 
synthesized SUMER Fe XIX line only observes plasma with temperature greater than 6 MK, which 
only exists in the heated loop. So we extract position S in panel (a) of Fig. [2] to compare with 
the black line in Fig. [3l which is located at the same position in the confined loop length. Panel 
(a) of Fig. 0] shows the temporal evolution of the synthesized AIA 131, 94 A channel emission, 
number density and temperature of position S, and panel (b) represents the temporal evolution 
of the synthesized Doppler shift velocity and intensity of SUMER Fe XIX line of the black line, 
which is close to the right footpoint. So the second and third peaks of light curve in panel (a) 
don’t mean peaks of the exactly second and third periods of waves, actually they indicate the time 
before and after the wave impact. The reason why the third peak is stronger than the second peak 
in AIA channels is that of AIA channels are remarkably affected by the density. Since after 
the impact, the front of wave-flow reflects back and collides with the tail of itself and the density 
at position S increases as shown in panel (a) as well as the third peak of AIA channels. We find 
that the number density, temperature and both AIA channels emission have an in-phase relation¬ 
ship as seen in panel (a), while the same is quantitatively true for Doppler shift and intensity in 
SUMER Fe XIX line in panel (b). Thi s suggests that these reflected patterns are propagating waves 
which show an in-phase relationship (jSakurai et al.ll2002l L rather than standing slow mode waves 
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which show a quarter-p eriod phase lag between velocity and intensity disturbances (|Wang et al 
2003al : lYuan et alJ 12015! ) . The difference between AIA 131, 94 A channel emission is because AIA 
131 A emission is more sensitive to higher temperatures around 10 MK, while AIA 94 A emission is 
sensitive to the temperature around 6.5 MK. This also explains the reason why AIA 94 A emission 
increases more than AIA 131 A emission at the third peak in panel (a) of Fig. O because at that 
moment the temperature decreases below 9 MK, while the density increases. 

In order to calculate the period of the reflected patterns, we extract another position C shown 
in panel (a) of Fig. [2j The position C is located at the apex of the loop, so the light curve of AIA 
131, 94 A emissions of slice C in panel (a) of Fig. [5] can reveal the more correct half period of the 
reflected patterns rather than panel (a) of Fig. 01 We identify three peaks of the reflected patterns 
based on the light curves of AIA 131, 94 A emission which all show strong damping afterwards. 
The time intervals between the three peaks are 285 and 360 seconds, so t he total period woul d 
be 570 seconds and 720 seconds. Compared with the observational result in iKumar et al.l (j2015l ). 
with 409 seconds for the period, the period in our simulation is longer. There are plenty of factors 
which could influence the period, such as the length and structure of the loop, the temperature 
inside the loop, etc. One possible reason for the longer period is the slower sound speed in our 
simulation, indicating that the reality should have even higher temperatures (greater than 10 MK). 
The period increase in our simulation is mainly affected by the large cooling, i.e th e strong thermal 
conduction which is mai nly responsible for the quick loop temperature decreases (lOfman &: Wane 


2002; 


Wang et al.ll2003bl l and little leakage of waves across the loop when the wave-flows impact the 


footpoint. When the temperature decreases during the propagation, and so does the sound speed, 
the period of the reflected patterns increases. The transition region at x = —23 Mm, shown by the 
AIA 304 A emission in panel (b) of Fig. [Sj also displays an oscillation. This oscillation indicates 
the height variation of the transition region, due to the impact from reflected patterns. The height 
variation actually traces the energy and mo mentum exchan ge between coronal and chromospheric 
regions as leaked by the reflected patterns. I Yu et al.l (|2013l l reports that quasi-periodic wiggles of 
microwave zebra pattern (ZP) structures can be associated with fast magnetoacoustic oscillations 
in a flaring active region. This oscillation of transition region associated with slow mode waves in 
our simulation may be observed in the future as well. 


4. Conclusion 


In this work, we performed a 2.5D MHD simulation to imitate the chromospheric evaporation 
and the following reflected patterns in a flare loop. We demonstrated that the periodic intensity 
variations ca ptured by the synthesized AIA 131 and 94 A emission images match well with previous 
observations (jKumar et al.ll2013l . l2015l l. 


With a particle tracer, we confirmed that these reflected patterns contain a clear wave com¬ 
ponent, in their sound speed like propagation. Through predicted paths of sound waves, we also 
found that these reflected patterns are dominated by the wave component while modulated by mass 

























- 10 - 


flows. To sum up, the reflected patterns observed in our simulation contain both slow waves and 
mass flows. 

With the synthesized Doppler shift velocity and intensity maps in SUMER Fe XIX line emis¬ 
sion, we confirmed that these reflected patterns are propagating slow mode waves rather than 
standing slow mode waves in our simulation, due to the in-phase relationship between Doppler 
shift and intensity. 

From the light curves of the synthesized AIA 131, 94 A emission, we estimated the period 
of oscillations which increases from 570 seconds to 720 seconds during the observed three periods. 
The increase of the period was due to the decreasing loop temperature and sound speed, caused 
by the strong cooling. The height variation of the transition region shown in the synthesized AIA 
304 A map may exhibit similar oscillations, correlated with the reflected patterns. This could be 
searched for in future observations. 

The research has been sponsored by an Odysseus grant of the FWO Vlaanderen. The results 
were obtained in the KU Leuven GOA project GOA/2015-014 and by the Interuniversity Attraction 
Poles Programme initiated by the Belgian Science Policy Office (lAP P7/08 CHARM). Part of the 
simulations used the infrastructure of the VSC - Flemish Supercomputer Center, funded by the 
Hercules Foundation and the Flemish Government - Department EWI. 
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Fig. 1.— Temporal evolution of number density (left column), temperature (middle column), 
and synthesized AIA 131 A emission (right column) images at t ~ 0, 83, 166 and 581 seconds, 
respectively. The arrows in the left column mark the velocities and directions of the local plasma. 
The white lines in panel (f) denotes a hxed loop, defined by a field line with a fixed width of 1200 
km. The cross in the panel (i) illustrates the initial location of the particle tracer ( x = 0 Mm, 
y = 25 Mm). An animation of this hgure is available. 
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Fig. 2.— Values of number density, temperature and the synthesized AIA 131 A emission inside 
the loop shown as time-distance maps are displayed in panel (a), (b) and (c) , respectively. The 
six red solid lines in panel (b) and (c) show the paths of virtual particles propagating at the local 
sound speed. The dotted line in panel (a) shows the temporal evolution of the actual Lagrangian 
particle tracer. Position S and C are used to analyse the light curves in Fig. [l]and[5l respectively, 
tl, t2 and t3 indicate the times at which the tracer particle changes its direction. There is another 
online animation of evolution of virtual particles in panel (a) of Fig. [2j 
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(a) Doppler Shift (100 km/s) 
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Fig. 3.— Doppler shift oscillations revealed by the synthesized SUMER Fe XIX line emission maps: 
Panel (a) shows time series of Doppler shift in Fe XIX along a top view slit and panel (b) shows 
time series of the Fe XIX line intensity. The black full lines show the location of position S used in 
Fig. [2] in this top slit view, and are used to analyse the light curves in Fig. 01 
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Fig. 4.— Panel (a) shows the temporal evolution of the synthesized AIA 131, 94 A channel emission, 
number density and temperature for position S in Fig. [2l Panel (b) represents temporal evolution 
of the synthesized Doppler shift velocity and intensity of SUMER Fe XIX line for the black line in 
Fig. El 
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Fig. 5.— Panel (a) shows the temporal evolution of the synthesized AIA 131, 94 A channel emission, 
number density and temperature for position C which locates at the loop apex as shown in Fig. [21 
Panel (b) a Time-Distance plot illustrates the height and intensity variation of the plasma at the 
left footpoint (x= -23 Mm) imaged by the AIA 304 Channel. 














